Scaling of V2 in heavy ion collisions 
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We interpret the scaling ol the corrected elliptic flow parameter w.r.t. the corrected multiplicity, 
observed to hold in heavy ion collisions for a wide variety of energies and system sizes. We use 
dimensional analysis and power-counting arguments to place constraints on the changes of initial 
conditions in systems with different center of mass energy y/s. Specifically, we show that a large 
class of changes in the (initial) equation of state, mean free path, and longitudinal geometry over the 
observed ^/s are likely to spoil the scaling in V2 observed experimentally. We therefore argue that 
the system produced at most Super Proton Synchrotron (SPS) and Relativistic Heavy Ion Collider 
(RHIC) energies is fundamentally the same as far as the soft and approximately thermalized degrees 
of freedom are considered. The "sQGP" (Strongly interacting Quark-GIuon Plasma) phase, if it is 
there, is therefore not exclusive to RHIC. We suggest, as a goal for further low-energy heavy ion 
experiments, to search for a "transition" y/s where the observed scaling breaks. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Nq 
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I. INTRODUCTION 

The azimuthal anisotropy of mean particle momen- 
tum (parametrized by it's second Fourier component V2), 
thought of as originating from the azimuthal anisotropy 
in collective flow ( "elliptic flow" ) , has long been regarded 
as an important observable in heavy ion collisions. The 
main reasons for this is that elliptic flow has long been 
understood to be "self-quenching" P, Q : The azimuthal 
pressure gradient extinguishes itself soon after the start 
of the hydrodynamic evolution, so the final V2 is insensi- 
tive to later stages of the fireball evolution and therefore 
allows us to probe the hottest, best thermalized, and pos- 
sibly deconfined phase. 

In addition, as has been shown in |3] , the V2 signature is 
highly sensitive to viscosity. The presence of even a small 
but non-negligible viscosity, therefore, can in principle be 
detected by a careful analysis of V2 data. 

Indeed, one of the most widely cited (in both the aca- 
demic and popular press) news coming out of the heavy 
ion community concerns the discovery, at RHIC, of a 
erfect fluid" , also sometimes referred to as "sQGP" 
[1, S 13 ■ The evidence for this claim comes from the 
successful modeling of RHIC V2 by boost-invariant hy- 
drodynamics [13i, [3, [iHl . The scaling of V2 according to 
the number of constituent quarks further suggests that 
the flow we are seeing is partonic, rather than hadronic, 
in origin [l, [§] , especially since the scaling applied to ki- 
netic energy (rather than transverse momentum) holds 
for every known species up to the lowest momentum 

While hydrodynamics is a fully deterministic theory, 
it contains a crucial not very well understood assump- 
tion: initial conditions. While the degree of boost in- 
variance is not currently well known experimentally, the 
transverse structure of the energy density should follow 
a Glauber model [l^, based on the superposition of the 
initial nuclear densities. That allows us to characterize 



the collision in terms of a number less than one called 
the eccentricity e, related to the impact parameter h and 
the radius R (~ A^/^fm) 
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the total transverse area of the system S depends, simi- 
larly, on R and b 



S = 2i?2 cos^i 



b 

2R 



b\ 



fe2 

T 



Since the announcement of the discovery of the perfect 
fluid, a considerable amount of high quality experimental 
data has been collected. In particular, extension of RHIC 
beams to smaller colliding systems such as Cu — Cu have 
allowed us to compare systems of similar multiplicity but 
in very different energy regimes. The results have been 
remarkable: It seems that W2/e (where e is the initial ec- 
centricity), plotted against (where S is the area of 

the collision system and ^ is the multiplicity rapidity 
density), fall on a "universal" curve, which links very dif- 
ferent regimes, ranging from Alternating Gradient Syn- 
chrotron (AGS) to RHIC ([i3,[il,Fig. [T|). 

This scaling, albeit with a linear dependence rather 
than the slightly curved one observed, has been predicted 
previously 19, 20] on the basis of a nearly "free stream- 
ing" calculation where the mean free path is compara- 
ble to the system size. Such a limit considerably under- 
predicts the observed W2/e, which is why a nearly-perfect 
hydrodynamic regime is thought to apply. As we will 
show, however, the same scaling in this regime is far from 
guaranteed. 

We do not possess at present the tools, such as 3D vis- 
cous hydrodynamics, to quantitatively analyze this data. 
However, the extent of the scaling suggests that these 
very different systems vary somehow only in one scale, 
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FIG. 1: (color online)A compendium of evidence for the uni- 
versal (for the observed energies) scaling of vije vs [l^ . 
Predictions from ideal boost-invariant hydro are also shown 
(Lines, see references of [T^). 



and that this scale is connected to the total entropy pro- 
duced [m. 

In this work, we shall perform a qualitative analysis, 
using elementary tools such as Taylor expansion, dimen- 
sional analysis, and power counting. We show that these 
tools, together with experimental data, allow us to place 
stringent limits on the initial conditions of the system 
created in heavy ion collisions, from AGS to RHIC ener- 
gies. 

In particular, we show that the observed scaling places 
very strong constraints on initial microscopic properties 
(entropy density, mean free path), as well as longitudi- 
nal structure. We argue that statistical and transport 
properties can not significantly vary between RHIC and 
observed lower energies. We therefore conclude that the 
"perfect fluid" , if it is there, is a common characteristic 
of all experimentally studied systems so far. 

It should be underlined that the "experimental result" 
on which we base our conclusions is, itself, somewhat 
theory-laden since e is not an experimentally measured 
quantity. If the best physical description of the soft initial 
dynamics at RHIC is not the Glauber model but rather, 
for example, the color glass condensate [13], the scaling 
might need to be revisited [l^ . Even within the Glauber 
model calculation, it is only when eccentricity fluctua- 
tions ^24J are taken into account that the universal scaling 
is observed (one of the reasons why this scaling was not 
reported until recently, and missed for example in ^25'|). 
Yet, the fact that an improvement in the computation of 
e results in improvement of the scaling is indicative that 
something physical must be behind it. 
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FIG. 2: (color online) Top panel (a): The range in 
XjjS^I'^ for the available experimental data Bottom panel (b): 
The range of S and e covered by available data 



It should be noted that the result in Fig. [5] is based on 
integrated vi, and is therefore subject to systematic 
uncertainties due to the wide range of acceptance in the 
experiments summarized in Fig. [2l Recent investigations 
(Fig. 1 of [2^, [13]), however, seems to show that the 
scaling is not dependent on integrating over since 
a compatible scaling holds when separate pt bins are 
considered. The universal scaling was found to hold in 
rapidity space also [13 

In the rest of the paper, we will regard the result as es- 
tablished because we wish to show that its consequences 
are profound. Finding simple scaling in a system as com- 
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plicated as a heavy ion colhsions is encouraging enough 
that it deserves exploration even through the scaling ob- 
servation should be regarded as preliminary. We hope 
that the experimental community will soon determine 
how universal this simple scaling really is. 



II. HOW NATURAL IS THE OBSERVED 
SCALING WITHIN HYDRODYNAMICS? 

For the system under consideration, three dimension- 
less parameters can safely be assumed to be significantly 
smaller than unity: The initial spatial eccentricity e, the 
initial speed of sound Cg (that encodes information on 
the equation of state [H, [2^, [s^l), and the initial mean 
free path Imfp divided by the initial size of the system, in 
general parametrized by the transverse radius R and the 
longitudinal size (z). The latter refers to the system's 
longitudinal size in configuration space, at mid-rapidity, 
at the start of the hydrodynamic evolution, and is a gen- 
eral definition. It is usually thought that (z) can inter- 
polate from a Landau initial condition [31|], where (z) is 
related to the nucleon mass m^r and the center of mass 
energy y/s by 



(z) ^ Ry^ 

rriN 



(3) 



, to the Bjorken [32!| initial condition, where, for a system 
with maximum rapidity y^, 
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where tq is the thermalization timescale of the system. 

Since hydrodynamic evolution is fully determined by 
the initial conditions and the equation of state^ , any flow 
variable can be thought of as a function of the parameters 
characterizing these. If this function is integrable (i.e., 
if no turbulence occurs), than it is safe to expand this 
function around any dimensionless parameter less than 
one. 

Both the dimensionless V2 and the dimensionful 
fm~^) can therefore be Taylor-expanded around 
these quantities. For V2, we know that the 0th term is 
(perfectly central collisions have no V2), so 
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An additional variable experimental observables can depend on 
is a freeze-out criterion, encoded by, for example, fr^i^i-ott ^ 



We disregard this variable in the subsequent discussion as the 
observables we discuss should not strongly depend on it. 
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aijk are in general (probably transcendental) functions 
of an arbitrary number of dimensionless quantities con- 
structed out of (z), T,iiB, R (in general ( = J2mni Cmn!, 
where Cm«; = (z)"* T"/i^i?'+"-™, aU m,n,l). The exact 
form of aijk can be obtained integrating the hydrody- 
namic equation from the initial time to freeze-out time 
(through they are expected to be insensitive to the latter, 
and equivalently to Il^iii^^y Ideal boost-invariant hy- 
drodynamics, with a bag model or lattice-inspired equa- 
tion of state m, [l3,[ll], predicts that 020^ ^ oioj, and 
hence W2/e is approximately constant (as the lines in Fig. 
[T] show) . 

The experimentally observed rise of W2/e with multi- 
plicity, and encounter with the hydrodynamic calculation 
(see Fig. [Ij , can therefore be interpreted as RHIC energy 
being the only point where the system reaches the "ideal 
hydrodynamics limit" . 

It is however unclear to what extent is such a conclu- 
sion an artifact of the models being used to perform this 
comparison assuming exact (isl . [l5| or approximate 
[33l . [s^ l boost-invariance as an initial condition at all en- 
ergies. While there are good physical arguments for why 
such an initial condition is appropriate for heavy ion colli- 
sions at mid-rapidity 32], the fact that some experimen- 
tal data is more compatible with Landau hydrodynamics 
even at RHIC highest energies [slj suggests the need to 
question this assumption, and in particular to evaluate 
it's effect on our estimate of the transport properties and 
their energy dependence. While no viscous calculation 
using Landau initial conditions has so far appeared in 
the literature, it is reasonable to suppose that the slower 
cooling from a Landau initial condition would leave more 
time for V2 to form. Hence, the limits on viscosity /mean 
free path/thcrmalization time inferred from V2 data are 
strongly correlated with the degree of assumed boost in- 
variance. 

contains information about both the longitudinal 
structure at freeze-out and the final particle number den- 
sity (a function of the initial T and /is). In an ideal 
(isenthropic) expansion, the final entropy is equal to the 
initial entropy content of the system (~ the initial parti- 
cle density n(T, /is)), so 



ldN_ 

S dy 



{z)n{T, hb) 



(7) 



Collective evolution, if the system has a non-negligible 
mean free path, can however create additional entropy. 

The first correction to isenthropic expansion should 
therefore be proportional to the entropy creation due to 
viscous processes. This is given by [36| 
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where 77 is the viscosity, Ui, the flow field, and V{t) the 
volume of the fluid. 
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viscosity is in turn proportional to the mean free path, 
density n, and mean momentum current (p) jsij 

V ~ lvafpn{T, ^b) (p) {T, hb) (9) 

since the initial volume is, to 0th order in eccentricity, 
~ (z) i?^,we get 

r-^ /^\„(T III,] + h'^ V/P , l{2>Wp, 

S dy -^(^)"i^'MBj(^l + 6oio ^ +^010 + 

Cs (6?n^+4^i^)+-) (10) 

where, once again, the coefficients b are functions of 
any arbitrary sum J2m n i (-z)™ T'"a*b-R"^'~'", that have 
to be calculated by integrating the viscous hydrody- 
namic equations. Since entropy is predominantly pro- 
duced in the initial collisions ^STJ , we can again disregard 

It is immediately clear that several of the parame- 
ters used in the previous expansions (in particular the 
Cs,lmfp, (z)) can not, by causality, depend on the trans- 
verse system size A, and have to depend only on the 
local energy density only. If soft observables scale with 
the number of participants, therefore, these parameter's 
dependence on energy and system size can only be a func- 
tion of a, where 

a = ^ (11) 

Not all parameters, however, have this dependence: The 
initial transverse system size R and eccentricity e exhibit 
no energy dependence. Rapditiy [27| provides a further 
indipendent direction. 

Presently available heavy ion experiments have ex- 
plored a significant range in e,cr and R, shown in Fig. 
m Had V2/e scaled in a different way from w.r.t. 
these variables, the currently available experimental data 
should have signaled it by the appearance of "branches" , 
systems with same but different V2/e (or vice- versa). 
Hence, the absence of such branches is a probe capable 
of constraining the nature of the system created in heavy 
ion collisions. In the rest of the paper, we will use this 
probe, as well as Eqs [6] and [TOl to constrain the system 
as much as possible. 

Fig. [3] sketches the dependence of Cs , Imfp and (z) 
as conventionally thought. As the energy increases, Cs 
rapidly dips in the "mixed phase" ^ and than increases as 
the system moves from a Hadron Gas to a Quark Gluon 
Plasma. An increase in Cs, if other variables were fixed, 
would increase the V2/e (since there is more pressure 



^ In a first order phase transition, it dips to zero. In a cross-over, 
it merely decreases j30| 



build-up) and, since it's an equilibrium process, either 
maintain constant or decrease it (if longitudinal flow 
is significant). 
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FIG. 3: (color onIine)The expected dependence, with energy 
of the initial (z) ,Cs and Imfp- These quantities, by definition, 
will be independent of system size 

Imfp is sensitive to the phase structure in a some- 
what different manner: the system should first enter the 
strongly coupled low viscosity sQGP regime, and then, 
when a is very high (y^ ~ TeV), the higher viscos- 
ity high- Temperature regime where the QGP becomes 
asymptotically free. A higher viscosity means a lower V2 
3, but a higher as more entropy is created in the 
system and collective energy is transformed, by micro- 
scopic interactions, into thermal energy. 

Finally, the stopping power of the system should de- 
crease as the initial condition goes from Landau [SH to 
the Bjorken [s^l limit. This decreases the pressure build- 
up needed to create V2 as well as -j^f^, since entropy is 
re-distributed in a wider rapidity space. The combined 
effect of these three scalings is difficult to evaluate with- 
out a solver of 3D viscous hydrodynamic equations. It is 
however likely to be non-negligible, and depend only on 
a, not on A or centrality/e. The total size of the sys- 
tem, and it's lifetime, however, should depends strongly 
on both a and A, as well as the eccentricity. 

Hence, terms ^ r Tr''{z) ' present in both Eq. [6] and 
1101 should vary in a non-trivial way as energy and sys- 
tem size are changed. In addition, for systems with a 
short lifetime, or a slow formation of V2, terms ~ e*^^ in 
Eq. [6] should become non-negligible. Thus, the appear- 
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FIG. 4: (color online) Scaling expected with for "typical" 
values of aijk,bijk and the formula in [ssl ] for multiplicity. 
The bands, lower to upper, correspond to different classes in 
centrality (10) and different system sizes (Pb,Au,Cu,S). 



The temperature is found by solving the resulting con- 
servation of energy equations, using the equation of state 
given in Eq. 1121 The resulting temperature is then 
plugged into Eqs. [14] and [13] to calculate Cg and Imfp- 
These, together with "typical" coefficients aijk,bijk and 
R ^ A^/^, (z) ^ To sinh(yi) are then used to investigate 
how U2/e and depend on each other. 

Since all "small" dimensionless parameters are encoded 
in Cs, e, ^0^, all coefficients a^jk ~ V2\^^p , hjk ~ 1. Fig. 
H] shows what kind of scaling is to be expected from Eqs 
[S] and [TU] if the approximate equalities hold exactly. The 
branch structure is clearly seen across the experimentally 
studied system sizes. It is of course possible to eliminate 
and produce universal scaling, but, in the absence of a 
deeper principle why that should be so, the coefficients 
a, b would need to be carefully fine-tuned. 

Thus, the experimentally observed scaling of Fig. [T] 
places very profound constraints on how the microscopic 
properties, and the global longitudinal structure, can 
vary between AGS and RHIC energies. 



III. WHAT DRIVES THE SCALING? 



ance of a common scaling between systems with different 
(7, A, e,such as the one observed appears highly unlikely 
a priori. 

We shall further explore the naturalness of the ob- 
served scaling using natural values of aijk,bijk, together 
with parametrizations of and Imfp 



p = ^(l + tanh(r,-r)) + 
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PQGP 



(l + tanh(r-rc)) 



(12) 



0.1 



Cs = — (l + 0.9tanh(T,-r)) 



Adeal 



(l + tanh(T-Tc)) 



(13) 



^ = ^(l + tanh(T,-T))+ (14) 
fm 2 

H (1 + tanh(r - rj) * log (^1 + J 

We then use the Bjorken formula to get the initial longi- 
tudinal and energy density distribution 



IdN 1 dN 



S dy ttA^ dy 



Top(T, hb) 



(15) 



where ^ can be obtained from experiment by a phe- 
nomenological formula jsl] . 

dN Npdi^iif.'ipants 



dy 
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The point made in the previous section is actually 
straight-forward to understand: As illustrated in Fig. \E\ 
a universal scaling means that the two quantities that 
scale are functions of a common variable (that can be, 
in general, a function of still further variables). 

Thus, the systems from AGS to RHIC appear to be 
controlled by a common scale, related to the total multi- 
plicity, which varies smoothly and drives both V2/e and 
■5 ■ This conclusion is a strong indication that micro- 
scopic properties of the system (equation of state and 
mean free path) are unchanged, up to a shift related to 
this scale, in the experimentally accessed energy range. 
It also raises the question of the exact nature of the vari- 
able that drives the scaling. 

It has been suggested, in [H,!!^, that the system, even 
at it's initial stage, is not entirely in the nearly inviscid 
QGP phase, but a fraction of it is in a highly viscous 
hadron gas. In this case, Eq. [6l and [TOl should be updated 
with new parameters Cs ^ ^QGP,hg ^ Imfp 
and an additional parameter 



iQGP,HG 
mfp 



(dN/dy) 



QGP 



(dN/dy) 



(17) 



total 



should be added. 

Could it be that a is what moves from zero to unity 
in the curve of Fig. [T]? If we assume a Glauber model 
and a Saxon- Wood distribution for initial density, and 
a critical "transition" energy density pc (independent of 
energy and system size), a becomes 

^ ^ / pjx, y, z, A, y/s)e{p{x, y, z, A, y/s) ~ pc)d^x 
J p{x,y,z,A,y^)d^x 
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and TA^sif) is the usual transverse participant density 
of the target nuclei. 



(20) 



Such scaling is very non-trivial to model exactly, but it 
seems to obey an approximate error function dependence 



(21) 



it is difficult to sec how such a scaling can be compatible 
with the observed universal curve. In particular, there 
will be limiting values in at which a should go at 
and 1 independently of A and h. above/below these 
scales, the system reduces to a one-component system 
and branches should start appearing. 

Similarly, at low energies there will be a regime where 
a = below a certain impact parameter, and a > above 
this impact parameter. In the first class of centralities, 
the large mean free path in the "corona" would introduce 
a dependence on A that does not conform to the univer- 
sal curve. While RHIC experiments have isolated several 
low centrality bins where vi is significantly below ideal 
hydro predictions (the corona dominates), such a criti- 
cal centrality has so far not been observed, (although 
this can be explored further with collisions of very heavy 
nuclei, such as Uranium [iol. fiTjV 

We conclude that, if the scaling parameter is indeed 
the QGP fraction of the system, the energies explored so 
far are far from the regime where this fraction is either 
close to unity or to zero. While this is possible, given the 
large variation in energy within the systems explored so 
to date, it would be somewhat surprising. 

An alternative ansatz to reproduce the observed scal- 
ing is by postulating a universality in initial conditions, 
up to a scale parameter, that also smoothly controls the 
initial temperature (and hence Cg and Imfp)- The basic 
constraint that the universal scaling imposes on a set of 
uniform initial conditions is the requirement that a sin- 
gle dimcnsionful scale (r) exists, and varying either the 
energy or the system size only shifts the system up and 
down the scale (r). Both (z) n{T, ^b) and V2 are then 
functions of only (r) as well as constants independent of 
energy and system size. A natural interpretation for (r) 
is the system's lifetime in the co-moving frame. In units 
of the mean free path, this corresponds to the inverse of 
the Knudsen number [425, the number of collisions be- 
tween the system's degrees of freedom. 



Because of the leading dependence of 



1 dN 
S dy 



(z) niT^fiB) (Eq. fTU]) . it follows that (r) can only scale 
with A, ^/s, a in the same way as (z) n(r, /is). 



{T)=F~\i,z)n{T,iXB)) 



(22) 



^ where F{{t)) is the same for all energies. A different de- 
pendence would lead to two different (r)s corresponding 
to the same (z) n{T, hb), which would again break the 
observed scaling. 

Eq. [S] and [TO] then simplify to 

V2 ^ e (aioo + aioiCs + aio2C^ + ...) + 

^mfp 



-e (flllO + fllllCs 



e (ai20 + ai2iCs + •■•) + 



(23) 
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&011 



s^mfp 



(24) 



where lmfp,b,a depend only on (r), presumably via a sim- 
ple scaling between (r) and initial temperature. This 
scaling will most probably be monotonic 



= F [ ^ 



T, 



ijiitial 



(25) 



More complicated scalings, with minima and sharp tran- 
sitions, will in general lead to a violation of the universal 
scaling, since events with similar final multiplicities could 
in this case have different W2/e. 

To fully appreciate these constraints, it must be re- 
membered that the lifetime strongly depends on the sys- 
tem's longitudinal initial conditions, and in particular 
Landau and Bjorken type initial conditions, with the 
same equation of state and transport coefficients, will 
lead to very different (t)s [43| . 

Thus, the scaling of V2 rules out a transition of the sys- 
tem from the Landau to the Bjorken limit such as the one 
seen in the top panel of Fig. [31 in the considered range of 
energies and system sizes: Such a transition would mean 
that two events with the same ^ ^ , one high-energy non- 
central, the other low-energy central, would correspond 
to two different lifetimes (the first close to the Bjorken 
limit, the second to the Landau limit), and hence, in gen- 
eral, to two W2/e. 

The monotonic increase of V2/e further constrains ei- 
ther the initial conditions to be far away from the Bjorken 
limit, or the mean free path to be non-negligible, at all 
energies : As initially inferred in [|| , and explicitly shown 
in p^.[T^.[T5j V2 is a self- quenching signature, which satu- 
rates after a finite time Ti,2, with r^2 <C (t) in the Bjorken 
limit [l^ 0, [il] ■ The V2 scaling than implies that the 
system never reaches r^2- If it did, systems with different 
{t)^2 ''"i'2 but the same e would have the same W2/e. 
It is clear that such systems would, in general, have very 
different -j^y^, breaking the scaling. 

The universal scaling of V2 in pseudo-rapidity space 
observed in adds a further layer of constraints to 
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Eq. [231It appears that the cr is also a uniquely de- 
termined universal function of the relative position of 
the volume element in rapidity space. This rules out 
a "Landau" /fire-strak model where the system is closely 
localized in rapidity. Instead, the fireball evolves in a 
way that is both local in rapidity, and strongly rapidity- 
dependent. Perhaps the BGK initial condition [46], could 
provide such an ansatz, although it would imply that such 
a geometry holds, to a good approximation, up to AGS 
energies. 

The picture suggested by the "universal" scaling of V2 
is then considerably different from the "RHIC reaches the 
perfect fluid" scenario: The "fluid" produced in heavy 
ion collisions, at all energies and system sizes where the 
scaling holds, should have a comparable Imfp, equation 
of state (cg) and longitudinal structure of the initial con- 
dition. All that varies is lifetime (r), which is uniquely 
determined by the initial density multiplied by longitudi- 
nal size, (z) n{T, ^b)- This universal structure is robust 
inasmuch the scaling experimentally observed. 

The alternative is that "we all have got it wrong" 
and the picture quantitatively analyzed in (loj , of a very 
weakly interacting system, is the appropriate one for de- 
scribing heavy ion collisions, from AGS to RHIC. But, 
aside from the difficulty in modeling such a large V2 in 
this picture, the conclusions in the previous paragraph 
would actually not be changed: For the scaling shown 
in [l^ to hold, it is necessary for {vijcrij), where Vij is 
the relative speed and cr^ the cross-sectional area of the 
system's microscopic degrees of freedom, not to change 
significantly with energy and system size. This is equiv- 
alent to requiring the microscopic properties of the sys- 
tem, such as viscosity and equation of state, to remain 
the same. 

Experimentally it will be very interesting to see at 
what point, in low energy collisions, is the universal 
scaling observed here broken. This point could well be 
the critical a that produces a deconfined system. 

Perhaps a greater energy and system size exploration 
around the region of the so called "kink" ," horn" ," step" 
[44,] anomaly can yield discoveries; As seen in Fig. [5] 
(lower panel, ~ 6/m~^), there might be a hint of 

splitting in the scaling curve; While the error bars abun- 
dantly drown out any firm evidence at this point, the 



approximate coincidence with the features highlighted in 
[4J| , as well as the breaking of HBT radii scaling with 
multiplicity and -^s [isj are suggestive. Likewise, it will 
be interesting to see if peripheral LHC collisions, where 
the initial temperature should lie solidly in the asymp- 
totic freedom regime, can be related to central top energy 
RHIC collisions. 

The observation that the scaling could apply not just 
with integrated pr but within pr bins ([2^, [23]) also 
deserves further exploration. One could object that, 
since {pt) is energy dependent, a meaningful compar- 
ison across different energy regimes can not be ob- 
tained. However, in the hydrodynamic picture (pt) 
is also strongly system-size dependent because of the 
growth of transverse flow, yet [la, [23] finds a pt specific 
scaling to hold when Cu-Cu and Au-Au are compared. 
We await further, energy dependent results in this direc- 
tion, but remark that the breaking of this scaling could 
signal the energy scale at which equilibration stops ap- 
plying. 



IV. CONCLUSION 

In conclusion, we have shown that the scaling of W2/e 
with places tight constraints on the hydrodynamic 
initial conditions in heavy ion collisions. It imposes an 
energy-independent relationship between initial energy 
density and longitudinal size, and makes it likely that 
longitudinal structure and microscopic parameters, such 
as the initial temperature, equation of state and viscos- 
ity are comparable in the considered range of energies 
and system sizes. We have suggested that looking for 
when the given scaling breaks might yield information 
about the critical energy and system size at which we 
can speak of a deconfined collective phase. 
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